Quantum tunneling of vortices in MgB2 superconductor 
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Abstract 

Magnetic relaxation in a MgB2 superconductor was measured. The temperature dependence of 
the normalized relaxation rate was determined for three different magnetic helds. By extrapolating 
these rates to T = K, we find that these extrapolations do not approach zero, indicating quantum 
tunneling of vortices in the superconductor. A quantum correction of the relaxation rate, followed 
by the correction of the magnetic moment, is proposed. Using the quantum correction, we find 
that Uq increases with decreasing temperature and approaches a maximum at T = K. 

PACS numbers: 74.60. Ge, 74.60.Ec, 74.40.+k, 74.70.Ad 
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Recent experimental studies of vortex dynamics of MgB 2 superconductor suggest that 
quantum tunneling of the vortices may strongly influence the superconducting behav- 
iors Thus, better undersdanding the mechanism is very important for future appli- 
cations. Magnetic relaxation technique is often employed to characterize superconduc- 
tors. During the measurement, a magnetic field is first applied to a superconductor at 
a fixed temperature, and then the magnitude of the magnetic moment m(t) is measured 
vs. time. The m(t) is frequently analyzed with Anderson-Kim (AK) flux-creep model ||, 
m(t) ~ ?7io[l — (T/Uq) ln(t/t )], where Uq is the energy barrier of the vortex motion in the 
absence of a driving force, to the reference time, and mo the magnetic moment at t = to- 
Deviations of the AK model are found in high T c superconductors (HTSCs), leading to a 
more general relation || |J 

m(t) » ^ (1) 

[l + ^T/Uo) ln(t/t )]i; 

where the characteristic factor /i = — 1 for the AK model, fi — > for the Zeldov model ||. 
and < /i for the collective creep model []3], [|, This relation predicts that the normalized 
relaxation rate 

_ dlnm T T ^ uT , , . S1 , . 

i2 = « — — w — 1 - — Ht t ) , (2) 

dint U + fj,T\n(t/t ) U [ U V 7 ° ;J ' V; 

[for the diamagnetic moment or for the persistent current (m oc j, where j is the current 
density)] vanishes linearly with decreasing temperature. However, experiments showed that 
the relaxation rates of superconductors do not extrapolate to zero ||, 0, |9], |Kj at T = 



K, indicating the decay rate being probably governed by quantum tunneling of vortices 

In this study, magnetic relaxation is measured at fixed temperatures in three different 
magnetic fields. We find that the extrapolations of the magnetic relaxation rates do not 
approach zero when T — > 0, suggesting the existence of quantum effects in the vortex system. 
Based on the theory of quantum tunneling vortices and the magnetic relaxation data, a 
quantum correction relation is proposed, followed by corresponding barrier evaluations. 

A small MgB 2 thin chip (thickness about 0.2 mm, and weight about 2.2 mg) was studied 
0| . The sample shows a sharp superconducting transition temperature around 39 K with 



a corresponding nearly single phase of MgB2 as determined by X-ray diffraction. Using 
a magnetic property measurement system (SQUID, Quantum Design MPMS 5.5 T), the 



FIG. 1: Magnetic relaxation data at magnetic field [XqH = 2.0 T for 10 different temperatures. 

magnetic relaxation of the sample was measured in magnetic fields with a scanning length 
3.0 cm. The measurements are processed at several fixed temperatures and magnetic fields 
at 1.0, 2.0, and 3.0 T. The sample was first zero- field cooled to a measurement temperature 
from a temperature (T = 50.0 K) well above the superconducting transition temperature 
T c , and then a magnetic field was applied for the measurement. 

Figure [l] presents magnetic relaxation data of the sample in the log-log scale at 2.0 T 
for 10 fixed temperatures from 2.0 to 22 K. All curves in Fig. |l] approximately linear in 
the log-log scale, except for the data in the initial time stage as illustrated in Fig. |l], which 
could mainly correspond to two reasons: the linear simulation (in the log-log window) being 
too simple for simulating each m(t), and reconfiguration of the flux line gradient, caused 
by a suddenly applied magnetic field. For a simple approach as shown with linear fits 
in Fig 1, we however can omit the data in an initial stage, and simulate the data using 
Zeldov model, mm M =° lnm — ln[l + (fiT/Uo) ln(t/t )]//i »° him — (T/Uq) ln(t/t ), with 
mo = m(to) and to — 1 s. The facility of this simulation is that leads to a simple relation 
R = —dlnm/dlnt ~ T/Uq without changing the qualitative conclusion of the energy barrier 
from a more precise simulation. 

Figure |2] presents the R(T, H) curves of the MgB 2 sample at different magnetic fields 
for t > 600 s (for largely reducing the influence of the reconfiguration of the flux line 
gradient in the initial stage). These data are strongly magnetic field dependent and show 
approximately linear relations in the linear-log window, resulting in that each R curve is 
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FIG. 2: Temperature dependent relaxation rates at magnetic fields 1.0, 2.0 T and 3.0 T. 



approximately simulated by a monotonic relation R = R exp(aT) as shown with solid 
lines, where Rq is the quantum creep at T = K, and a the magnetic field dependent 
constant. Note that this monotonic relation gives a remarkable comparison to that of a 
HTSC for which a broad peak or a plateau of a R curve in the intermediate temperature 
range (0.2T C < T < 0.8T C ) is observed |g |, §. R{T,H)s of a HTSC show monotonic 
reductions with decreasing temperature at low temperatures, where quantum tunneling of 
the vortices may play an important role for the magnetic moment decay. Considering the 
similar R decreasing tendency over whole superconducting temperature regime for MgB 2 , 
we may conclude that quantum tunneling exists through the regime. 

According to the quantum tunneling theory || ITTJ, [l2| for the single vortex tunneling 
in the limit of strong dissipation, the decay rate R{T) « h/S(T) ps e 2 Pn{T) /%L c (T) = 
Q{T)\j c {T) / jo{T))^ 2 , where S is the saddle-point solution (bounce) of the Euclidean action, 
L c = C-Jjo/jc the collective pinning length along the magnetic field direction, Q = e 1 p N jfi^ 
the dimensionless quantum sheet resistance, p^ the normal state resistivity, £ the coherence 
length, jo the depairing current density, and j c the critical current density. As a consequence, 
we can assume that the relaxation rate is the sum of two contributions corresponding to 
the purely thermal contribution, and the additional quantum tunneling contribution respec- 
tively. Hence, the relaxation rate is: 



FIG. 3: Energy barriers determined by different processes. The solid circles, and open triangles 
correspond to the Uq determined by Uq = T/R, and Uq = T/(R — Rq) respectively. The solid 
lines, and dash lines represent the relations Uq = T/Rq exp(aT), and Uq = T/i?o[exp(aT) — 1], 
respectively. 

where m q is the magnetic moment using the quantum correction, R = R(T = 0) and 
f(T) accounts for the temperature dependence of the quantum contribution, relating to the 
temperature dependent S(T). Integrating Eq.(|3|), we obtain the magnetic moment for the 
quantum correction: 

m (t) » ^ (4) 

t*°f[l + (piT/U )]n(t/ta)}» 

Considering a typical experimental time t < 10 5 s and a typical decay rate Rq < 0.01 with 
< f(T) < 1, one will find the additional factor t Rof < 10 05 p» 1.1 leading to m(t) w m q (t) 
and d In m/d In t ~ d\nm q /dlnt in a broad temperature regime and a limited time scale. 

The interest is the comparison of the barrier evaluations between Eq.(fJ) and Eq.(^). 
The solid circles in Figure |3| show the U (T,H) data determined using Eq.(|2|) with \i = 
(Uq = T/R). The solid lines represent the corresponding fits of the data using the relation 
Uq = T/R = T/Rq exp(aT). From the fits, it is easy to find that each Uq(T, H) curve shows 
a maximum value (peak) in the intermediate temperature regime and Uq(T — >) —>■ in the 
very low temperature regime. The tendency of Uq(T — > 0) — » is unreasonable and contrasts 
sharply with the other superconducting properties, e.g. the lower the temperature, the larger 
the irreversibility field and the larger the critical current. For the barrier evaluations with 
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Eq.(^|), the temperature independent contribution, corresponding to f(T) = 1, is taken 
]T3, 0, leading to the barrier U = T/(R - R f) = T/(R - Ro) as shown in Fig. § with 
the open triangles. Remarkably, using the approach R = R exp(aT), we find that each 
Uq = T / RQ[exp(aT) — 1] curve (dash line) in Fig. monotonically increases with decreasing 
temperature and reaches its maximum value Uo = l/R^a at T = K, which notably 
contrasts to the corresponding Uo = T/R = T/R Q exp(aT) curve (Uq{T — > 0) — > 0). For the 
extrapolations, one will find that the smaller the Rq, the larger the quantum correction in 
low temperature; besides, one will find that quantum correction is relatively small in high 
temperature where thermal effects are dominant. As a result, we conclude that a quantum 
correction of the relaxation rate is necessary for properly evaluating the energy barrier. 

According to the quantum mechanics, the tunneling probability of vortices depends on 
the barrier height. Since the height is temperature dependent, the f(T) = 1 may not be in 
reality. Therefore, we assume that both R(T) and f(T) can be simulated by polynomials, 

R = Ro + £ R n T n (5) 

and 

/ = 1 + £ f n T n , (6) 

n>l 

where R n and /„ are polynomial coefficients, and n is the integer for both polynomials. 
Taking Eqs.(]5|) and (|5]) to Eq.(^), we have 

U (T) ^T(R-Rof)- 1 

w (E RnT"- 1 -RoY, fnT"- 1 )- 1 , (7) 



n>l n>l 



leading to 



In general, we have 



U o (T = 0)n(R 1 -R f 1 )- 1 . (8) 



(R-R )/T=J2RnT n - 1 >0. (9) 

n>l 

Using the boundary conditions: < f(T) < 1, f(T = T c ) = and f(T = 0) = 1 , we 
determine 

(/ - 1)/T = Y fnT^ 1 < 0. (10) 

ra>l 
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Using Eqs.@ and (|10D to Eq. (0), one will find that Uq(T) increases with decreasing 
temperature and reaches its maximum at T = 0. 

In summary, we have measured the temperature dependence of the magnetic relaxation 
rates of the MgB 2 superconductor in three magnetic fields. Based on the theory of quantum 
tunneling of vortices and the magnetic relaxation data in this work, a quantum correction 
of the relaxation rate, followed by the correction of the magnetic moment, is proposed for 
evaluating the energy barriers of the vortex motion. Using this correction, we find that Uq, 
which sharply contrasts to the result without using the correction, increases with decreasing 
temperature and approaches a maximum at T = K. 

This work has been financially supported by the National Science Foundation of China, 
and PAI 4/10 (Belgium). 
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